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The dynamics of a dense relativistic quantum fluid out of thermodynamic 



■ equiHbrium is studied in the framework of the <I>^ scalar field theory in the 

large N limit. The time evolution of a particle distribution in momentum 
space (the tsunami) is computed. The effective mass felt by the particles in 
such a high density medium equals the tree level mass plus the expectation 
value of the squared field. The case of negative tree level squared mass is 



Q_i' particularly interesting. In such case dynamical symmetry restoration as well 

^ • as dynamical symmetry breaking can happen. Furthermore, the symmetry 

D , may stay broken with vanishing asymptotic squared mass showing the presence 

' of out of equilibrium Goldstone bosons. We study these phenomena and 

identify the set of initial conditions that lead to each case. We compute the 
^ ■ equation of state which turns to depend on the initial state. Although the 

■ system does not thermalize, the equation of state for asymptotically broken 

symmetry is of radiation type. We compute the correlation functions at equal 
times. The two point correlator for late times is the sum of different terms. 
One stems from the initial particle distribution. Another term accounts for the 
out of equilibrium Goldstone bosons created by spinodal unstabilities when 
the symmetry is asymptotically broken. Both terms are of the order of the 
inverse of the coupling for distances where causal signals can connect the two 
points. The contribution of the out of equilibrium Goldstones exhibits scaling 
behaviour in a generalized sense. 
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I. INTRODUCTION 



Physical phenomena in high energy density situations cannot be treated with the usual 
perturbation methods. Self-consistent non-perturbative methods are necessary in order to de- 
scribe the out-of-equilibrium dynamics of relaxation of quantum fields in such situations. The 
large limit is a particularly powerful tool for scalar models. The need for a self-consistent 
method stems from the fact that the particle propagation in such situations depends on the 
detailed state of the system. 

Relevant physical situations at high energy densities arise in the hadron (quark-gluon) 
plasma as in ultrarelativistic heavy ion collisions, in the interior of dense stars and in the 
early universe 

We consider a dense relativistic quantum gas out of thermodynamic equilibrium. We 
investigate the non-perturbative dense regime in which the energy density is of the order 
m|j/A where A is the scalar self-coupling and m^j the physical mass in vacuum. In these 
conditions, even for very weakly coupled theories non-linear effects are important and must be 
treated non-perturbatively. We thus use the large N limit that provides a non-perturbative 
scheme which respects the internal as well as the space-time symmetries, is renormalizable 
and permits explicit calculations. More precisely, we consider the 0{N) vector model with 
quartic self- interaction and the scalar field in the vector representation of 0{N). 

We study the physics of the in-medium effects at high energy density. In such out of 
equilibrium regime the effective mass felt by the particles is time-dependent and different 
from the tree level mass m^. In the large- iV approximation, it takes the form 

M\t) = ml + ^{$'){t) (1.1) 

where stands for a A^-component scalar field with self-coupling A. 

We see that the quantum fluctuations of the scalar field directly contribute to the effective 
mass. We shall consider for simplicity translationally invariant situations. That is, invariant 
under spatial traslations. Therefore, the effective mass and all one-point expectation values 
depend on time but not on space coordinates. 

We work in the small coupling regime A <^ 1. The reason being that the different 
dynamical time scales are well separated in the A <^ 1 regime. For larger couplings the 
dynamical time scales become of the same order and different physical phenomena get mixed 
up. 

We choose as initial conditions a gaussian wave functional and ($)(0) = ($)(0) = 0. The 
calculation of correlation functions is done in the general case ($)(0) 7^ 0. 

As initial distribution of particles we will choose a 'tsunami' 0. That is, a distribution 
of particles in momentum space that we choose for simplicity spherical. Typically, such 
distribution is a shell peaked in a wavenumber ko with a large density of particles ~ mj^/ A. 

We have seen that these conditions do not determine completely the initial state, and 
that the remaining degrees of freedom can be interpreted as the quantum coherence between 



2 



different fc-modes. We choose two highly coherent initial states that we will call case I and 
case II. 



As we see in eq. ( p,.l[ ) the in-medium effects give a positive contribution to the effective 
mass. In general, this contribution to the effective mass initially oscillates with time. These 
oscillations can produce parametric resonance which is here shut-off after a few oscillations 
by the damping of the oscillations due to decoherence. This resonant effect is negligible for 
small A. (Parametric unstabilities are shut off by the backreaction for no particle initial 
conditions ||-|ig). 

We present analytic solutions for initial particle distributions narrow in momentum space. 
Such self-consistent solutions express in close form in terms of elliptic functions, and repro- 
duce the numerical solutions with very good accuracy till the damping of the oscillations 
becomes important (due to decoherence phenomena). The narrower is the initial distribution 
in momentum, the longer in time this effective solution holds. 

For m|j > 0, the effective mass is larger than the tree level mass > (see table 

1). 

The more interesting case corresponds to < 0. We have in such case spontaneously 
broken symmetry for low density and low energy states. For large initial energy density, the 
positive definite term |($^)(t) may overcome the negative squared mass m|. < in eq. ( |1 . 1| ) 
and the symmetry may be restored: A^^(t) > (see table 2). 

This happens at t = for the case we call II provided the energy density E fulfills 

For late times the symmetry is restored in both cases, I and II, provided 

|4 1^ 12 



^R ^R 

Moreover, we find in case II an interval of energies where the symmetry is initially unbroken 
(A^^(O) > 0) and where it is broken for asymptotic times: 



Aij Xr Xr 

That is, we have a dynamical symmetry breaking for these situations (see table 2). 

We compute the asymptotic equation of state for this out of equilibrium system. That is, 
we derive for late times explicit formulae for the pressure as a function of the total energy. 
The equation of state we obtain explicitly depends on the initial state, A and m^j (see table 
3). Notice that even for asymptotic times the system does not thermalize. In particular, the 
distribution functions reach non-thermal limits for t ^ oo. 

When the symmetry is asymptotically broken, the equation of state takes the radiation 
form 

P = ^. 
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in spite of the fact that the system is out of equihbrium. Moreover, the Goldstone theorem 
is vahd here in this out of equihbrium situation. Namely, the effective squared mass A4'^{t) 
asymptotically vanishes when the symmetry is asymptotically broken. 

The equation of state for asymptotically unbroken symmetry [eqs. ( [7.7|) ] has the cold 
matter and the radiation equation of state as limiting cases. 

We explicitly compute the two point correlation functions. 

For late times, the two point correlator at equal times t) expresses as a sum of two 

or three terms: 



C{\x\,t)=Cor^g^n{\x\)+Cp{\x\,t)+Cs{\xlt) (1.2) 

There is the time- independent piece Corigin{\x\) concentrated around the origin. The pulse 
term, Cp{\x\,t), is due to the particles in the initial distribution that effectively propagate 
as free particles with mass A4oo- 

Cp{\x\,t) = -l- F{\x\-2vgt-c) , 
X\x\ 

where Vg is the group velocity [vg < 1 for unbroken symmetry and Vg = 1 for dynamically 
broken symmetry where Aioo = 0]. F{u) is non-zero only around u = 0. 




The last term in eq. ( |1.2| ) is only present for dynamically broken symmetry. 

Cs{\x\,t)= ^ 

where i^' is a constant. The function Q{u) is of the order 0{1) only for < m < 1 due to 
causality. When the order parameter is identically zero we have Q{u) = 9{1 — u). 

When the order parameter is nonzero, Q{u) oscillates with u. At a given time t, the 
number of oscillations of Q{u) in the interval < m < 1 equals the number of oscillations 
performed by the order parameter ($) {t) from time t = till time t. That is, scaling exists for 
($)(0) 7^ in a generalized sense since the function Q{u) changes each time performs 
an oscillation. This is due to the appearance of an extra length scale, the initial value of the 
order parameter. 



II. THE MODEL 

We consider the 0(A^)-invariant scalar field model with quartic self- interaction with 
the scalar field in the vector representation of 0{N). 
The action and Lagrangian density are given by, 

S" = y d^x C , 
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The canonical momentum conjugate to is, 

n(x) = (2.2) 

and the Hamiltonian is given by, 

H = Jd'x |ln^(x) + 1 [V^x)r + lrn'^' + ^ («f } • (2-3) 

In the present case we will restrict ourselves to a translationally invariant situation, i.e. 

eigenstates of the total momentum operator. In this case the order parameter < t) > 
is independent of the spatial coordinates x and only depends on time. 

The Heisenberg equations of motion for the field operator take the form 

= (2.4) 

The coupling A is chosen to remain fixed in the large limit. 
It is convenient to write the field in the Schrodinger picture as 

$(x) = (a(x), ff{x)) = (ViV0(t) + x{x), V{x)) (2.5) 

with {f]{x,t)) = where ff represents the A^ — 1 'pions', and (j){t) = (o"(x)). Thus, (xi^)) = 0. 



A. The wave functional (Schrodinger picture) 

We shall consider Gaussian wave functionals of the type 



^^.,t)]=^^'/'it) n.^exp 



where 



Vki^) = d X r/(x, t) e 



ik-x 



(2.6) 



(2.7) 



[Hence, we can assume A_j:{t) = Aj^(t) without loss of generality]. 

As shown below, such Gaussian wave functionals are solutions of the Schrodinger equation 
in the N = oo limit. The Hamiltonian (|2.3|) in the large N limit is essentially a harmonic 
oscillator Hamiltonian with self-consistent, time-dependent frequencies: 



H{t) = NVK,it)-—{j:{%-U) 



{21 



+ E 

k 



where hci{t) stands for the classical Hamiltonian 



+ o ^K^) ■ v-k 



hciit) = - 



A 



and 



tu5(t) = P + + 



A 



The functional Schrodinger equation is then given by {h = 1) 



(2.9) 



(2.10) 



More explicitly, 



(2.11) 



which then leads to a set of differential equations for ^g(t)- 

The evolution equations for Aj^{t) and A/'(t) are obtained by taking the functional deriva- 
tives and comparing powers of on both sides. We obtain the following evolution equations 



Af{t) = Af{0) exp J f dt' 



2NVhd{t') 



AN 



k Ik 



(2.12) 
(2.13) 



with A^^{t) = A^^{t)^%A^^^{t). 

The time dependence of the normalization factor Mif) is completely determined by that 
of the A^{€) as a consequence of unitary time evolution. 

Using the expression for the wave functional (|2.6|) we have, 



'<nk-v-k)-- 
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iV 



2 A^dt) 



and 

leading to the self-consistency condition, 



N 



NV 



or taking the infinite volume limit, 

(r/2(x)) 1 r d^k 



N 



NJ (27r)3 



(2.14) 



(2.15) 



(2.16) 



The expectation value of the Heisenberg equation of motion 
motion for the order parameter 0(t) [Q: 



0(t) +m2 0(t) + ^ 



02 (t) 



N 



m = 



(2.17) 
yields the equation of 

(2.18) 



where we used that (ff) = and (f/^ ff) = 0. 

Eqs.( |2.12D -( p.l3| ) together with eq. ( |2.18| ) define the time evolution of this quantum state 
in the infinite N limit. 



B. The field modes (Heisenberg picture) 

The Ricatti equation ( |2.12| ) can be linearized by writing Aj:{t) in terms of the functions 

^j:{t) as 

-Mt) = , (2.19) 



leading to the mode equations 



n + Vl = 0. (2.20) 



The relation ( |2.19| ) defines the mode functions ifj^ up to an arbitrary multiplicative con- 
stant that we choose such that the wronskian takes the value. 
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2i . 



(2.21) 



The functions (fj^ have a very simple interpretation: they obey the Heisenberg equations 
of motion for the pion fields obtained from the Hamiltonian ( |2.3| ). Therefore, we can write 
the Heisenberg field operators as 



rf(f , t) 



v^(27r)= 



ik-x 



+ 4 e 



-ik-x 



(2.22) 



where , at are the time independent annihilation and creation operators with the usual 
canonical commutation relations. Thus, the fj^{t) are the mode functions of the field. 
From eq. ( |2.19| ) we obtain the following useful relations, 

1 . , . Id 



^Rkit) 



Ajkit) 



2dt 



(2.23) 



Then, using these relations and eq. ( p.l7| ) we express the fluctuation in terms of the mode 
functions: 



N 



d^k 
(2^ 



(2.24) 



C. Definition of the particle number 



Since in a time dependent situation the definition of the particle number is not unique, 
we choose to define the particle number with respect to the eigenstates of the instantaneous 
Hamiltonian ( |2.8|) at the initial time, i.e. 

1 52 



Nujj, 



2 5% ■ 6f]_^ 

Here, Uj^ is the frequency (p^ ) evaluated at t = 0. 

The expectation value of the number operator in the time evolved state is then 
where A^(t) and 5g(t) are defined through the relations. 



(2.25) 



\kit) 



(2.26) 



(2.27) 



In terms of the mode functions v^^(t) and ^%{t), the expectation value of the number 
operator is given by 



4:UJr 



(2.28) 



For initially broken symmetry the frequencies in eq. (p.25|) are modified for low k according 
toeq. ( pJTI) . 
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D. Initial conditions 



We will take initial conditions for the field expectation value of the form, 

0(0) = 00 , 0(0) = , 
since we can make 0(0) to vanish by a shift in time. 



(2.29) 



The initial quantum state defined by eq. ( pl6|) is determined giving y4^^(0) and ^/^(O) 
for all k or equivalently, using the relations ( |2.27| ), giving = A^(0) and 5g = 5^(0). 

^m(0)=^dl + ^^ . AikW=^k^k- (2-30) 
We can invert eq. ( ^.2(j| ) expressing Ag(t) in terms of %(t) and the particle number n^{t), 



(2.31) 



Thus, we can express the initial conditions in terms of the initial distribution of particles 
n^(0) and 5^. However, there are two possible values of A^ for each n^(0) due to the ± sign 
in eq. (|2.31|) . We shall call case I to the upper sign in eq. (|2.31|) and case II to the lower 
sign. (It is very important to notice that this sign is not the sign of m\. But is due to the 
fact that the value of 7^ and 5^ do not fix completely the initial state for the modes.) 

The quantity 5^(t) appears as the phase of the wave functional and will be chosen for 
simplicity to be zero at t = 0, 



(2.32) 



For simplicity, we will consider spherically symmetric particle distributions and with gaussian 
form: 



ng(0) = — exp 



{k - ko) 



(2.33) 



where Nq is the total number of particles per unit volume and / is a normalization factor. We 
shall always consider Nq ^ rn?. As we shall see below (subsection |IV A|) it will be convenient 
to consider Nq ~ m? / g. 

As the initial conditions are spherically symmetric and the evolution equations are in- 
variant under rotations the solutions will be spherically symmetric. So, the dependence on 
k is only through the modulus k. 



The initial conditions on the mode functions follow from the relation ( p.l9|) and the initial 
conditions ( |2.3m ) on A^{t) plus the wronskian constraint ( |2.21 ). Thus, the initial conditions 
on the mode functions are: 



</^fc(0) 



vjfc(O) 



(2.34) 
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with rik defined by: 



= ARk{0) = Uk [1 + Afe] 



(2.35) 



One sees that for 6k - 
conditions become: (pk{0) 
used in refs. iBHTDI . 



and an initial state with no particles (A^ = 0), the initial 
-- Xj ^Ju\ and 0fc(O) = —i^fuTk- These were the initial conditions 



III. EVOLUTION EQUATIONS 

The evolution equations for the expectation value of the field 0(t) and for the mode 
functions ^Pkit), eq. ( |2.18| ) and eq. ( p.20| ) respectively, are: 

^ + + ^ + iv')B{t)]^ m = , (3.1) 

+ [k' + ml + ^ + {v')s{t)]^ Mt) = , (3.2) 

with the self-consistent condition 

The above evolution equations must be renormalized. This is achieved by demanding 
that the equations of motion be finite. The divergent pieces are absorbed into a redefinition 
of the mass and coupling constant. 

Mlit) ^ml + ^ [{r^')s{t) + 0^(t)] = ml + ^ [(r^^)^(t) + <P'{t)] . (3.4) 

A detailed derivation of the renormalization prescriptions requires a WKB analysis of the 
mode functions tfk{t) that reveals their ultraviolet properties. Such an analysis has been 
performed in refs. In summary, quadratic and logarithmic divergences are absorbed 

in the mass term while the coupling constant absorbs a logarithmically divergent piece 
The renormalized quantum fluctuations take the form 



2 J (27r)3 

with K, an arbitrary renormalization scale. 



(3.5) 



We define now dimensionless quantities choosing the physical mass in vacuum \mR\ as 
unit of mass: 
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k 



9 = 



\mR\ 



\mR\ 



2\mR\ 



(3.6) 



In terms of these dimensionless quantities the equations of motion become 



|^±l + CV)+,S(r) 



C(r) = 



+ q'±l + C\T)+g^{T) 



1 , e{q-l)Ml{T) 



2q^ {turI'^ 



(3.7) 
(3.8) 
(3.9) 



where we have chosen the renormahzation scale k = \mR\ for simphcity and the sign ± 
corresponds here to the sign of m\. 
Ineqs. (^-(il), 



r = 



\mR\^ 



« + C'(r) + (7S(r) 



(3.10) 



plays the role of a time dependent effective squared mass, a = signim^^) = ±1. 

Depending on whether A^^(O) is positive or negative the symmetry will be initially un- 
broken or broken. 



The initial conditions become in dimensionless variables for A^^(O) > 0: 

C(o) = Co , C(o) = o, 



0,(0) 



(3.11) 



(3.12) 



with: 



nq=UJq[l+^q] , Uq = Jq^+M^iO) 



A, = 2 



n,(0)±^n2(0) + n,(0)-| 



(3.13) 



We call case I to the upper sign and case II to the lower sign. (Notice that this sign is 
not the sign of m^. We have to distinguish cases I and II due to the fact that the value of 
ng{0) 7^ and 6g do not fix completely the modes.) 
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From ( p:32|) and (g^) we have: 



Sr, = 



with A/'n = and a 



Therefore, 



(3.14) 



1 , 9{q-l] 



2g3 



MHt) 



(3.15) 



To perform the numerical evolution, we have introduced a momentum cutoff A. The 
quantum fluctuations become for finite cutoff. 



1 - flog A [Jo 



q^dq\^,iT)\' 
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A^ + ^WA 



O 



g_ 

A2 



. (3.16) 



This is a positive quantity [up to 0{gy\. a = signlmj^) = ±1. 



In the case where the symmetry is initially broken {Ai'^{0) < 0) the only change in eqs. 
( p.lll )-( ^l^ ) is that the initial frequencies Uq are modified for low q as follows f^: 



g2 + |A^2(0)| for g2 < 1 + \M\0)\ 
g2 + 7^2(0) for q^ > 1 + \M\0)\ 



(3.17) 



IV. EARLY TIME DYNAMICS 

We study now the time evolution of a narrow particle distribution nq{0) peaked at g = go- 
We show below that we can approximate the dynamics for early times using a single mode 
ifefflr). We solve the time evolution of <feff{T) in close form in terms of elliptic functions. 
Moreover, we compare these analytic results with the full numerical solution of eqs. ( ^.7[ )- 

(PD- 

For < 0, we only consider initial particle peaks with gg > 2. Thus, they are well 
outside possible spinodally resonant bands. 

For simplicity, we consider = in this section and in sections M, and [VI1|, 



A. Case I 



In this case we have: 



We see from eq. ( |3.13| ) that t 



nq{0) + ^nl{0)+nq{0) 
le modes with nq(0) -C 1 have Q„ 



ujq. Therefore, v^g(O) 

and 0g(O) are of order [see eq. ( p.l2|) ]. Thus, their contribution to g'^ij) will be of order 
g for early times. 
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On the other hand, we have for modes with ^^(0) = 0{l/g) (as g <^ 1 this imphes 
n,(0) > 1) 



Ag ~ 4 72^(0) > 1 



fi, ^ 4 cjg n,(0) > 1 . 



Therefore, eqs. (|3.12|) imply that these modes have |v2g(0)| <^ 1. 
Thus, 



^?S,(0) 
A^'(O) 



0{g). 
a + 0{g) 



±l + 0{g) 



(4.1) 



(4.2) 
(4.3) 



where a = sign(mjj) = ±1. 

For these modes with nq(0) = 0{l/g), |0g(O)| ^ 1. Therefore, these modes will then 
grow and their contribution will dominate (y'S(T) for early times. 

We always consider that ^^(0) = 0(1/ g) for some interval in q. Hence, its contribution 
to fl'S(r) will be of order one and will have an important effect on the dynamics. In such 
conditions, the total number of particles A^^o per unit volume is also of the order 0{l/g). 

For the initial conditions considered (|3.14| ), these dominant modes are in a small interval 
centered at g = go? and they are in phase at least for short times. More precisely, the g-modes 
will stay in phase for {q — q^) t 27i . The modes which are in phase contribute coherently 
to 5'S(r) (each one with a contribution proportional to its occupation number). 

Hence, for small r a good approximation is to consider that all the particles are in a 
single mode with q = qo- This approximation will apply as long as the modes with large 
occupation number stay in phase. 

In these conditions eq. ( p.l6| ) yields. 



9 ql l^effir)]^ Ag 



Oig) + Oiga) 



(4.4) 



where Aq 
Thus, 



a [see eq. ( P.14|) ] 



dr2 



+ ^go "Peffir) +g Aqql \^eff{T)\'^ ^eff{r) = 



(4.5) 



with UJ, 



90 



go + « [recall a = sign(m|.) = ±1 and we choose q^ > 2]. 



We also have, 



A^o = 47rgo Ag Ueff 



No CO, 



eff 



go 



nq^Aq 



(4.6) 



Solving the non-linear differential equation ( [4.5|) with the initial conditions ( |3.12 ) and 
BD, we obtain (peff^r). We find from eq. 



(7S,(r) = (go^ + a) 1 



1 



1 - 2A;2 



1 — sn^ 



gg+g , 



(4.7) 
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with sn(2, k) the Jacobi sine function and 



k 



\ 



,1 I- 



(4.8) 



The function ( [4 .71) is non-negative and oscillates between zero and 

max = (go + 

with period 




2gNo 
TT (g2 + a)3/2 



(4.9) 



T 



2./V^ ir(A:) 
V go + « 



(4.10) 



where i^(/c) stands for the complete elliptic integral of the first kind. [Notice that eq. 
implies 1 - > 0]. 



The elliptic solution ( [4.7|) correctly predicts the amplitude of the first oscillation and the 
oscillation period (if the initial distribution of particles is not too wide around g = go)- See 
Fig. I. 

The amplitude of the numerical solution is well reproduced by the elliptic solution ( |4.7| ) 
till damping becomes important. The oscillation period keeps well reproduced by eq. ( |4.7| ) 
for longer times than the amplitude. [See Fig. ^ . We see from the numerical solution of the 
exact equations eqs. (|3.7|) -(|0|) that 5'S(r) exhibit significatively damped oscillations after 
a few periods. 



The integral over q for 5'S(r) [eq. ( |3.15| )] gets damped with time due to the lost of 
coherence between the different g-modes in the distribution peak. We can apply here the 
adiabatic approximation (see Appendix In the late time limit, 5'S(r) therefore tends 
to the value g'Ejmax/'^- In addition, the narrower is the peak, the slower the oscillations in 
5'E(r) are damped. See Fig. |^. 

The time scale where the numerical and the early times solution deviate in amplitude, 
essentially depends on the width a of the initial distribution. The smaller is a, the latest the 
early time solution ( |4.7| ) holds. Notice that eq. ([4.7|) gives g^^r) to zero order in cr, whereas 
the damping of the oscillations is given by higher orders in a. 

In addition, we have seen from the numerical resolution that the early time evolution 
depends on g only through gNo as predicted by eqs. (|4.^ )-( ^^ . 



When there is particle production through parametric resonance, particles in the initial 
peak distribution are annihilated (in order to conserve the total energy). This reduces the 
contribution of the initial peak to 5'S(r) whereas the oscillations due to the created particles 
give a contribution to (7E(r). These oscillations are due to the coherence between the created 
particles at different g. 
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Such changes in the distribution of particles influence the asymptotic value of 5'S(r) for 
late times. Contrary to the vacuum initial conditions, parametric resonance shuts off here 
by the damping of the oscillations and not due to backreaction as in ref. Therefore, 
for small g the influence of parametric resonances is highly suppressed. When parametric 
resonance is appreciable, it makes the dynamics to depend on g and not only through the 
combination gNQ. 

In case I, as Ai'^{0) = a = sign(m|,) = ±1. The symmetry is initially spontaneously 
broken or not depending on whether m\ (the squared tree level mass) is negative or positive. 



B. Case II 



In this case: Ag = 2 



As in case I, 5'S(r) is dominated for short times by the modes with large occupation 
numbers which are in phase at small r due to the initial conditions ( p.l4|) . Thus, we can do 



the same approximation as in Case I, considering that all particles are in a single mode with 
Q = 

We have for the modes with q qq, ng{0) = 0{l/g) ^ 1, 

A„ ~ -1 + , fi„ ~ < 1 . (4.11) 

^ 4n,(0) ' 4n,(0) ^ ' 

Thus, for r = 0, 

gHim = ^ +0{g) + 0{ga), (4.12) 

M\<d) = a + gi:u{^) , (4.13) 

where a = sign(m|.) = ±1. 

Eq.( |4.12| ) is a third degree equation in gT,ii{0) since ujg^ = \JqQ + a + 5'S//(0). The 
explicit solution is given in Appendix D. We find in the limiting cases 



gNo ) 



gEniO) ' "^-^^^^ ' (^^j - -(go' + a) + 
gEuiO) ^--<i-)^^^ + O ( -P^y (4.14) 



For early times eq. (|3.16|) becomes 

gE{r)=gqlAq\^,jf{T)\^ +0{g) + 0{ga). (4.15) 

Therefore, 
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+ {qI + «) ^effir) +gql Aq |</?e//(r)P y?e//( 



Using the initial conditions given by ( p.l2|) and 5e//(0) = we get from ( [4.11|) 



1^, 



where we have used that rie/f = 4^^°^^ - 



(4.16) 



(4.17) 



Thus, the solution of eq. ( [4.16|) with the specified initial conditions can be written as 
gEn{T)=gJ:j{T + T/2) 



(4.18) 



(Qo + a) 1 + 



1-2/^2 



sn^ 



l-2fc2 



(r + T/2), k 



where sn(2;, k) is the Jacobi sine function, T is the real period 



T = 2J\^ K{k), 
V ?o + a 

K{k) stands for the complete elliptic integral of the first kind and 



k 



\ 



1 



1 I 9^11 max 

H „;! I „ 



(4.19) 



95- 



The expression for gT^nmax is here 

g'^iimax = g^ii{0) 



9N0 



TrJqi + a + gEiiiO) 



(4.20) 



Notice that the relation gT^nij) = gT,i{T + T/2) is true for given values of go, « and k. 

Here 5'S(r) oscillates between zero and gT^jjmax- While the effective squared mass A4^{t) 
oscillates between its initial value /^"^{O) = ±1 + gT,jjmax and its tree level value ±1 at the 
minima of 5fS(r). 

We see that this approximation gives us correctly the amplitude of the first oscillation 
and the oscillation period; but not the damping of the oscillation that is due to the dephasing 
of the initial particle distribution. The broader is the initial particle distribution, the more 
effective the dephasing works and the faster the damping occurs. 

As in case I, the smaller is a (width of the initial particle distribution), the later the early 
time solution ( [4.18| ) holds (see figures 2 and 3). Recall that the damping of the oscillations 
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is given by higher orders in a while eq. ( |4.18| ) gives gT^ij) to zero order in a. These higher 
orders in a will also break the relation: gT,ji{T) = gT,j{T + T/2); that we have find at zeroth 
order. 

In case II, A^^(O) = a + gT^nmax- [Recall a = sign(m|j) = ±1.] As gT^umax > 0, for 
> the symmetry is initially unbroken. Instead, for mj^ < 0, the symmetry is initially 
unbroken for gT^umax > 1- For gT^umax < 1 the symmetry is initially spontaneously broken. 



V. INTERMEDIATE AND LATE TIME DYNAMICS 

We discuss here the intermediate and late time behavior of the quantum fiuctuations and 
the effective squared mass. We mean by late time, times later than the spinodal time and 
than the damping time. We consider (q = 0. 

We observe from the numerical results the following common features in the late time 
behavior for a wide range of initial conditions: 

The asymptotic constant values of the magnitudes (mass, pressure, number of particles) 
depend on g, for small g, only through the combination gNQ. (Except when parametric 
resonance is important). 

Energy is conserved to one part in 10^ confirming the precision of our numerical calcula- 
tions. 

We have also seen that the mass tends to its limiting value oscillating with an amplitude 
that decays at least as ~ 1/r. A similar asymptotic decay has been found in this model for 
zero particle initial conditions and C(0) 7^ H]. 



In section VII we derive the asymptotic equation of state. 



A. m| > (a = +1) 

We have for cases I and II that, 

M\T) = l+gJ:{T), (5.1) 

although g^ir) has a different expression in each case. 
As 5'S(r) > the symmetry is unbroken for all r. 

The adiabatic approximation holds (see Appendix ^ since parametric resonance is neg- 
ligible in the weak coupling regime considered here. Hence, the asymptotic value of g'^ir) 
is and the asymptotic squared mass goes to (see Fig. P 

Ml = l + ^-^>0. (5.2) 
This result is in good agreement with the numerical calculations. 
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We find that the mass tends to this value oscillating with an amplitude that decays at 
least as ~ 1/r. 

The initial peak of particles becomes lower and wider, and the total number of particles 
slightly decreases compared to its initial value. 



m| > 


Case I 


Case II 




+1 


+ 1 + gTjjirnax 


initial 

symmetry 


unbroken 


unbroken 




1 ^ 1 9^1 max 


1 2^ 1 9^11 max 


late time 
symmetry 


unbroken 


unbroken 



TABLE 1. Initial and late effective mass and symmetry for m|j > 0. 



B. m| < (a = -1) 

We have for both case I and case II that, 

M'{T) = -l+gJ:iT). (5.3) 

Recall that in case I the symmetry is initially broken since 2/(0) = 0. In case II, ^//(O) = 
^11 max &nd the symmetry is initially broken (unbroken) for g^iimax ^ 1 {g^iimax 

We have two different asymptotic regimes: 

• 3^]^ax ^ ^. Asymptotically unbroken symmetry. 

In this regime the results are similar to those in the previous subsection V A , with 

Ml = -l + ^-^>0. (5.4) 
This is indeed the value we obtain numerically (see Fig. 
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We find that for small g there is no appreciable spinodal resonance here. This is so 
because A^^(r) oscillates around the positive value A^^, although for some intervals 
of time M'^{t) < 0. 

In this regime the symmetry is restored for late times due to the presence of a high 
density of particles. 

In particular, we have already noticed that for case I, A4'^{0) = —1. This does not 
change appreciably the dynamics and the symmetry gets restored provided sIl^s^ > x_ 

• s^HiiHL < I Asymptotically broken symmetry. 

In this regime for intermediate times (times earlier than the spinodal time r^) A^^(r) 
oscillates around the negative value, 

= -1 + < . (5.5) 

giving rise to spinodal resonances. (See figs. |^, |^ and ^) The dynamics turns to 
be similar as for a constant squared mass — /i^. That is, the spinodally resonant band 
is in the g-interval from q = to q = fi and the spinodal time Tg is the same as for a 
constant squared mass of — yU^ (see appendix 0). 



We can further distinguish: 

— when gT^max < 1; A^^(t) is always negative for times r < r^. 

— when gT,max > 1, A1^(t) can be temporarily positive due to oscillations. This 
happens in case I where 7W^(0) = —1. 

In case II, A^^(O) = -1 + gT^max > 0, but after a time of order T/2 , M'^{t) < 0, 
and A^^(r) continues to oscillate around the negative value given by eq. (|5.5|) . 
Thus, in case II we have dynamical symmetry breaking. 



At time the spinodal resonance has created enough particles (of the order 0(1/ g)) 
to give an important contribution to gTi^r). This finally makes (7S(r) to oscillate 
around 1. Thus, A1^(r) oscillates around zero and the spinodal resonance stops. The 
particles created by the spinodal resonance are coherent. This gives new oscillations 
to A1^(r). These oscillations get damped and the squared mass goes for late times to 
its asymptotic value, 

Ml = 0. (5.6) 



The vanishing of the effective mass is accompanied by the presence of Goldstone bosons 
out of equilibrium as in 
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Particles created by spinodal resonances remain in the g-interval from to /x. This 
creation of particles depletes the initial peak of particles keeping the total energy 
conserved. 

For T < Tg the dynamic depends on g, for small g, only through the combination gNQ. 
However, depends explicitly on g. 



m\ < 


Case I 


Case II 


7^2(0) 


-1 


— 1 + gT^IImax 


initial 

symmetry 


broken 


if 9^ 1 1 max < 1 

broken 


if 9^ 1 1 max > 1 

unbroken 




if alii^ < 1 



if alii^ > 1 

1 1 O"^! max 

2 


j-p 9^11 max ^ 2 




j-p 9^11 max \ 2 

1 1 9^11 max 
2 


late time 
symmetry 


broken 


unbroken 


broken 


unbroken 


change of 
symmetry 


no change 


dynamical 
symmetry 
restoration 


no change 


dynamical 
symmetry 
breaking 


no change 



TABLE 2. Initial and late effective mass and symmetry for mj^ < 0. 



VI. ENERGY 



After renormalization, taking dimensionless variables and introducing a momentum cut- 
off, the renormalized energy can be written as H 



2\mR 



4 



Eren = C, (6.1) 

e = -C'(r) + -C'(r) + -C'(r) + i^ 
2sv; 4^^^ 4 2 
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+ f/%'rf?[l0.Mr + a^J(r)|^,(r)r] -i [^S(r)]^ + 0{g) 



(6.2) 



where a = sign(m^). 

One can easily check that the energy is conserved using the renormahzed equations of 
motion (0)-(|3: 



Using the initial conditions (|3.11| ), ( p.l2|) , ( |3.14|) and Co = we obtain for the energy at 
r = 



1 1 — a g 
e = - 



4 2 



(6,3) 



This equation gives the energy in terms of the initial data. 
1 1^ 

4 2 



e — i is always positive because we consider Co = and Co = 0. 



A. Case I 



Let us consider an initial narrow distribution of particles which is peaked at g = go- 
Thus, we can use eqs. ( [4.1|) and (|6.3|) and consider all particles in a single mode with q = go 
as in previous sections. (Recall that in this case = + a.) 



The energy is then given by: 



i% + ") + 7 



1 1 - a 



2tt 



4 2 



0{g) + 0{ga). 



(6.4) 



B. Case II 



Under the same approximations [now using eq. ( [4.11|) ], the energy is given by: 



e={ql + a) 



+ \^-^ +0{g) + Oiga) 



(6.5) 



where cu^^ = q^ + a + 5'S(0) and gJ^u^O) is given by the relation ([4.12| ), i.e., gT,ii{0) 
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C. gTimax in terms of the energy 



The energy and gT^max have different expressions in cases I and II [see eq. ( [4.9| ) vs. eq. 
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However, it is important to notice that gT^^ax has the same expression in terms of the 
energy both for cases I and II, 



1 1 - a 



e 



4 



(go + a) 



(6.6) 



We have verified that eqs. ( |6.4|) -( |0|) are valid numerically for initial particle distributions 
with width a < 1. 



VII. EQUATION OF STATE 

We derive here the equation of state {i.e. the pressure as a function of the energy) for 
asymptotic times. As shown below, the asymptotic equation of state depends on the initial 
state. 

We consider Co = and 6q = 0. Notice that Co = and Co = implies C(t) = for all r. 



A. Sum rule 



A further expression for the energy follows by evaluating the right-hand side of ( |6.1|) in 
the r — i> cx) limit. Using eqs. ( [BlUD and ( |B13| ) yields 



1 1 - a 



+ j\UqM'^{^)+Mlg^o.-\{g^oo)' +0{g) 



(7.1) 



[The cosine terms in eqs. ( [B10| ) and ( |B13| ) for late time are fastly oscillant in q and thus 
they do not contribute to the g-integral in the infinite time limit.] 

Equation ( [7. 1| ) allows us to express the integral over (oo) in terms of known quantities: 
the initial data and A^^, 



\' dqM' [oo] 



/o - - ^ ' ■ 4 2 

This sum rule holds for Co = and 6q = 0. 



e--^-MlgJ:^ + \igi:^r +Oig). 



(7.2) 



B. Pressure and equation of state 

The renormalized pressure can be written as 
2\mn\^ 



Pren (^) 



p(r) 



A 



R 



p{t) 



-^ + C +9 q dq 
Jo 



0{g) 



(7.3) 
(7.4) 



22 



We analogously evaluate the pressure in the infinite time limit with the result, 



4 

Poo = -e + - dq M^{oo) + Ml, gT^^o • 



Using now the sum rule ( [7 .21) and = a + gT^oo yields 

1 1 1 — a 

Poo= 3e--^7Loo--^-. 



(7.5) 



(7.6) 



[Recall a = sign(m|,).] 

We consider narrow distributions of particles centered at g = go- To obtain the equation 
of state, we make the approximation of considering all particles in a single mode with q = q^, 
and we express the pressure as a function of the energy. 

The equation of state is the same for case I and case II. 

We have to distinguish between unbroken and broken symmetry for r = oo: 

• Ml>0 

In this regime (7E00 = ^^^""^ with gYj^^ax given by eq. ( |6.6|) . Thus, the equation of state 
( [7. 61 ) becomes 



Pc 



1 a 
3 ^~ 6 



'(g2 + a)2 + 4 e- 



1 1 



a 



(go + a) 



1 1 - a 
3 2 



(7.7) 



This equation gives p^o as a function of e for an initial momentum distribution centered 
at go. Therefore, the equation of state explicitly depends on the initial conditions. 

Let us consider some limiting cases. 

— In the limit, e — ^ <ti (go + a)^, we have the equation of state. 



1 / 1 1 - «\ / a \ 1 1-a 



For a = +1 (m|j > 0) this equation reduces to 

and interpolates between a cold matter (for go <^ 1), and a radiation (for go ^ 1) 
equation of state. 

In the opposite limit e — | ^ (gg + a)^ > 1 we have a radiation type equation 
of state Poo = I 
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This can only happen if a = — 1 (m|j < 0). Ai 
and we have a radiation type equation of state, 



2 

oo 



—1 + gT,ao = imphes (7S00 = 1 



Pc 



(7.10) 



The pressure is continuous at the boundary between broken and unbroken symmetry, but 



its derivative with respect to the energy have a discontinuity | -^r^- [See eq. ( [7.10| ) vs. eq. 



We see that the equation of state is determined by the energy, the momentum go around 
which the initial peak of particles is centered and the sign of the physical mass in vacuum 
m|j. But it does not depend whether we are in case I or II. 



Ml>0 



m\ < 



m| > 



equation 
of state 



Pc 



Pc 



+ 



<7o-l 



1 I 



Pc 



+ 



</o+l 



1 + 



4e 



- 1 



TABLE 3. Asymptotic equation of state in the different situations 



VIII. CORRELATION FUNCTIONS AND BOSE CONDENSATE 



The equal time correlation function is given for an arbitrary time r by 
{^\r,r)r,\^,T)) - (r/'^(f,r)) (r^^(0,r)) = 6^^' C{f,T) = 6^'' J |^,(r)p e^^^-\ (8.1) 



The initial conditions are specified in eqs. ( |3.11|) - (|3.14|) ; and as they are rotationally 
invariant, 



1 f°° 

C{r,T) = — qdq |<^9(r)p sin(gr) . 
In this section we consider Co = and Co 7^ 1 with Co <^ 1- 



(8.2) 
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A. Early time 

We have the same results for Co = and for Co ^ 1- 



1. Case I 

At r = the correlation functions is of order one, because |v5g(0)p ^ 1. [ |y9g(0)p ~ 1 
for the non-occupied modes and |v5g(0)p <^ 1 for the highly occupied modes, see eq. (^TTI).] 
For the highly occupied modes we have |0g(O)p = 0{l/g). Thus, for early times r = 0{1) 
these modes will have |(y9g(r)p = 0{l/g). This makes the correlation function be of order 
1/ g near the origin for early times, as we see in Fig. ^ 

2. Case II 

In this case at r = 0, the modes have |v5g(0)p = 0{l/g). This makes the correlation 
function be of order 1/ g near the origin for r = 0, as we see in Fig. ^ and 0. 

B. Late time 

The late time behavior of the correlation function depends whether the symmetry is 
dynamically broken or not. However, its behaviour is the same for both cases, I and II. 
We have to distinguish the two regimes: 

1. A^^>0 

In this regime we have the same results for Co = and for Co 1- 

We observe that at intermediate times a spherical pulse develops. This spherical pulse 
propagates with a constant radial speed given by the radial group velocity for q = qo and an 
amplitude that decreases as 1/r. The radial width of the pulse, L, remains approximately 
constant (see Figs. and ||). 

The group velocity is asymptotically given by. 



diOqoo 

^9 - 



dq 



q=<io ^ 



(8.3) 



We recall that here Ml^ = a + ^^f^ > 1 [see eqs. (U) and {^^]. [a = sign(m|) = ±1.] 
Asymptotically, the correlation function becomes the sum of two terms 

C(r, r) = Coriginir) + Cp{r, r) , (8.4) 
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where Corigin{i^) is the correlation function near the origin. This term is asymptotically 
t ime-indep endent . 

The pulse contribution to the correlation function, Cp{r, r), has approximately the asymp- 
totic form, 

Cp{r,T) = ^P{r-2vgT-c) , (8.5) 

here c is a constant of order one, and P{u) is of the order 0{1) only for —L/2 < u < L/2 
where L is the width of the pulse {i.e. the pulse is localized around r ~ 2t>gr + c). 

The pulse term is due to the particles in the initial distribution that effectively propagate 
as free particles. This is so since the effective mass in the mode equations (|3.2|) becomes 
asymptotically constant and hence the modes effectively decouple from each other. The 



pulse term is absent when there are no particles in the initial state |T0[ . 

In summary, the correlator is of order 0{1/ g) ^ 1 for 2vg t+c—L/2 < r < 2vg t+c+L/2 
whereas causality makes it to fall to values forr > 2vgT + c + L/2. 

2. Ml, = ^ 

In this regime the symmetry is broken and there were spinodal resonances for earlier 
times. 

Since the effective mass vanishes asymptotically the mode with g = behaves as 

<^o(r) =™ L^Kt (8.6) 



The Wronskian guarantees that neither of the complex coefficients L, K can vanish ||10|| . 
This linear growth with time can be interpreted as an out of equilibrium novel form of Bose 
Einstein condensation. We analyse below the contribution of this condensate Cgir^ r) to the 
correlation function. 

We have studied both Co = and Co 7^ (with Co ^ !)• 

a. Co = 0. For late time (r r^), the particles created due to the spinodal resonance 
contribute to the correlation function with a term Cs(r, r) of order \jg when r is in the 
interval (0, 2r). This term decays as 1/r. 

There is in addition a pulse term, Cp{r,T), in the correlator that moves away from the 
origin with unit velocity (remember eq. (|8.3| ) and that = 0). See Figs. [^-|TT. 

Thus, the correlation function is asymptotically given by 

C{r, r) = Coriginir) + Cs{r, r) + Cp{r, r) (8.7) 

Here, Corigini^), the correlation near the origin is asymptotically time- independent. 
The contribution of the pulse, Cp{r, r), has the form, 



1 

gr 



Cp{r,T) = — P(r -2t - c) 
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The contribution from the particles created by spinodal resonance, Cs(r, r), has the form: 

C,{r,T) = — Q(pi . (8.9) 



gr \2t 

where i^' is a constant and Q{u) = 6{1 — u). 

Introducing the dynamical correlation length ^{t) ~ 2r and defining the variable u = 
, Cs{r, t) can be written in the form, 

C.(r,r) = — ^^(l-n). (8.10) 
Using the customary notation for the scaling regime 

^^(^-^)= [^(^)|2(i-.) -^H- (8-11) 

with the anomalous dynamical exponent z = 1/2 (the naive scaling length dimension of the 
field is 1). In this case the scaling function is given by, 

I{u) = — e{l-u). (8.12) 
g u 

An analogous spinodal term Cs(r, r) in the correlator have been obtained for initially 



broken symmetry and no particles in the initial state in |10|. 

The spinodal term Cs(r, r) can be interpreted as follows. For times r later than there 
is a zero momentum condensate formed by Goldstone bosons travelling at the speed of light 
and back-to-back. That is, massless particles emitted from the points (0, r) and (r, r) form 
propagating fronts which at time r are at a distance t — Tg from the origin and from r, 
respectively. These space-time points are causally connected for 2(r — r^) > r. Otherwise, 
Csir, t) is not of order 1/g but of order one. 

An alternative interpretation of the causality step function goes as follows. Signals are 
emitted at the speed of light from all points in the condensate. We have causal connection 
between the points and r once signals starting from a given point arrive to both points. 
The earlier this happens is for the signals emitted from the half- away point at r/2. These 
signals need a time r = r/2 to reach both points. Hence, the correlator is of the order 1/g 
for 2 r > r. 



Moreover, the analytical derivation of eq. ( p.l2| ) from the low-g behavior of the mode 



functions given in ref. [|T^ also applies here. Notice that the scaling contribution of the 
Goldstone bosons to the correlator is different from a free massless propagator. We have 
here a 1/r falloff whereas a free massless scalar field has a 1/r^ falloff. 
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b- Co 7^ (with (o <^ 1). For late time (r ^ r^), the particles created by spinodal 
resonance give a contribution Cs{r, r) of order 1/g to the correlation function in the interval 



< r < 2r. See Fies. KMB 



The correlation function is asymptotically given by 

C(r, r) = Coriginir) + Csir, r) + Cp(r, r) , (8.13) 

where the time-independent piece Corigin{r) is the correlation near the origin. Cp{r, r) is the 
contribution of the pulse, 

Cp(r,r) = — P(r-2r-c), (8.14) 
gr 

and Cs{r, r) is the contribution from the particles created by spinodal resonance, 

^^(-'^) = 1^(7)^2(1^ (8-15) 

with ^(r) ~ 2r the correlation length, u = and z = 1/2. 

Ir{u) is no longer given by eq. ( ^.12| ) when Co 7^ 0. Qt{u) = ulriu) now oscillates with 
u. At a given time r, the number of oscillations of Qt{u) in the interval < m < 1 equals the 
number of oscillations performed by the order parameter ({t) from time r = till time r. 
That is, the scaling exists for Co 7^ in a generalized sense since the function Qt{u) changes 
each time C(''") performs an oscillation. This is due to the appearance of the extra length 
scale Co- See Figs. 13-16. 

As for unbroken symmetry, the pulse term is due to the particles in the initial distribution 
that effectively propagate as free particles. 

Both for Co = and Co 7^ causality makes the correlator of order 0{1) for r > 2 r + c. 
For r < 2t + c the correlations are of order 0{l/g) ^ 1. 



APPENDIX A: ADIABATIC APPROXIMATION FOR THE MODES 



In this Appendix we use the adiabatic form of the modes to evaluate the quantum fluc- 
tuations S(r). 

The modes '^qir) can be represented as M 



1 



-i dx Vq{x) _|_ ^ g« /o^ dx Vq{x) 



(Al) 



where and hq are constants and Vqir) depends on time. Inserting eq. ( |A1|) into eq. (|3.8| ) 
yields the following non-linear differential equation for Vqi^r) 



W 3 fVqiry 



2 Vqir) 4\Vq{T) 



+ V'{T) = q' + M\r) 



(A2) 
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The initial conditions ( |3.12| ) combined with the Wronskian conservation yields 



1 . 



As long as Vgi^r) is a real function, we have for the squared modulus 



flgl + \bq\ + 2\aq bq\ COS [2 J dx Vq{x) + 



(A3) 



(A4) 



where aq is a time independent phase. 

When q belongs to a parametric resonant band Vq{T) gets an imaginary part. We consider 
here the case where such unstabilities have a negligible effect. 

We are interested in the modes with large amplitudes \aq\ ^ 1 , ^ 1. That is, those 
in an interval of width a around the peak momentum go- For such modes, thanks to eq. 
we can approximate \aq\ ~ \bq\. 



Inserting eq. ([A4|) into the integral ( |3.15|) for S(r), and approximating the slowly varying 
factor by its average on a period (that we will denote i ) yields 



2 / IsP 



1 + COS ( 2 / dx Vq{x) + aq 



(A5) 



It follows from here the bounds 

< S(r) < 4 



For late times the integral of the oscillating cosinus in eq. ( |A5| ) vanishes. Therefore, 

(f dq . ,2 



S oo 



and 



S(CX)) — ^ ^max ■ 



(A6) 



APPENDIX B: SLOWLY VARYING PARAMETERS IN THE LATE TIME MODE 

FUNCTIONS 



Let us define the following slowly varying parameters 

1 



Mr) 



2 
1 



goo 



Bq{r) ^ - e 



<^g(r) + 



qoo 



(Bl) 
(B2) 
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with. 



uj^^^Jq2 + Ml, ; Ml^M\^) . (B3) 



These slowly varying parameters are asymptotically constant. 
We can express the mode functions in terms of them as follows: 

(^g(r) = Ag(r) e^"'-" + 5g(r) e"^"'-^ . (B4) 

Thus, the squared modulus of the modes is, 

|^,(r)|2 = |A,(r)|2 + |fi^(r)|2 + 2\A^{t)B,{t)\cos[2u,^t + </),(r)] , (B5) 



where 



The constancy of the wronskian implies: 

|2 I A 



A,(r)i?,(r)* = |A,(r)5,(r)|e^^^M. (B6) 



B,{r)\'-\A,{r)\' = ^, (B7) 



^qco 



plus terms with derivatives of the slowly varying parameters that vanish asymptotically. 
We define a slowly varying modulus [0] , 



M,{T)^^g^\A,{r)\^ + \B,{r)\^ . (B8) 

The main contributions to the physical quantities comes from the modes with occupation 
numbers of order for these modes |74g(r)| and |i?g(r)| are of order Therefore, 
Mg{T) is of order 1 for these modes. Moreover, eq. (|B7| ) implies that, 

mr)\' = \A,{t)\' [1 + 0{g)] , (B9) 

and we can thus approximate the squared modulus as follows: 

gMT)\^ = M,(r)2 {1 + cos[2u;,oor + 0,(r)]} [1 + 0{g)] . (BIO) 

We obtain from ( p^) a similar formula for the squared modulus of the derivative of the 
mode function, 

+ig(r) 6^"^'°°^ + Sg(r) 6-^"^'°°^ . (Bll) 

Using Aq{oo) = Bq{oo) = and the same procedure we have used to derive the equation 
for the squared modulus, we obtain: 

g\^qiT)\' = uol^ Mli^) {1 - cos[2c.,oo r + 0,(r)]} [1 + 0{g)][l + 0{r-^)] . (B12) 

Or equivalently. 
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Mr)\' = {gWq{r)\' - 2 Ml{r) cos[2uOq^ r + 0,(r)]} [1 + 0{g)][l + 0{r-')] . (B13) 
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APPENDIX C: CALCULATION OF THE SPINODAL TIME IN THE PRESENCE 

OF THE TSUNAMI 



Here, we obtain the early time solution for the modes in the spinodally resonant band, 
and we estimate the spinodal time, r^, both for cases, I and II. 

Recall that there is spinodal resonance provided mfj < and gT^max < 2. 

Let us assume that the time scale for the development of spinodal unstabilities, r^, is 
longer than the time scale of damping of the oscillations in 5fS(r), r^. This happens when 
the particle peak is wide enough. 

Before entering on the calculation of r^, let us remark that the numerical calculations 
show that the main effect of 5'S(r) for r < (and small g) is giving a constant positive 
term S^hls^ tj^e effective squared mass. In fact, this turns to be true in general, not only 
when the oscillations of g'^ir) are damped before r^. 

We have for Td < r < Tg-. 

-^^^M,^,/r)^l + ^<0. (CI) 
Hence, an approximate equation for the modes reads, 

(|^ + g^-/^^)^.(r)=0. (C2) 

Thus, the modes with q in the interval between and /i are spinodally resonant. 

We consider initial peaks of particles centered well outside the possible resonant bands. 
Thus, for a sufficiently narrow initial particle peak there are no particles in the spinodally 
resonant band. Therefore the initial condition for modes in the resonant band are: 



= 77^ = + l-^'(O)l)''^' , 0.(0) = = -^ {q' + |A^'(0)|)'^' . (C3) 



The solution of eq. ( |U1D for these modes is: 

, . 1 



• (C4) 



+ (V^' - q' + t,Jq^ + \M\0)\j e-^v^ 
We thus obtain for the squared modulus neglecting the exponentially decreasing term 

,2 



4(/i2-g2) /g2+|y^2(0)| 



4/i2 (l - ^) ^q2 + IM 



1 + /^^ e^-V-^^. (C5) 
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The contribution of the spinodal band to 5'S(r) is given by, 



S,(r)= / q'dq\v,{r)\\ (C6) 
Jo 

Inserting eq. ( |C5|) into eq. ( |C6| ) we obtain an estimation for the spinodal growth of the 
quantum fluctuations. 

To approximately evaluate this integral, we can make further simplifications in eq. (|C5|) . 
As < 1 and the contribution of the modes with q ^ fj, is exponentially suppressed, we can 
expand in to second order in the exponential and to zeroth order in the factor outside 
the exponential. 

|<^,(r) |2 ^ e^-^ e-"^ . (C7) 



In addition the integrand has its maximum at g = 0(0. and < /i < 1. Therefore we 

can do the approximations 1 + yU^ ~ 1 and \J q^ + |7V1^(0)| ~ 1 (for both case I and case II). 
Thus, 

|^,(r)p^^e2^'^e-^^ . (C8) 



Then the integral over q takes the value 



The spinodal time is by definition, the time where the unstabilities are shut off for 
all < g < yU. This happens when the spinodal modes contribution ^^(r) compensates the 
initial (negative) value of M^jj-{t) [see eq. (|C1|) ]. 

(7E,(r,)^/i2 (CIO) 

Therefore Tg is given by the following implicit equation, 

16 



+ ^log(/ir,). (Cll) 
4/i 



The spinodal times given by this equation are in good agreement with the numerical 
results. 



APPENDIX D: THE INITIAL QUANTUM FLUCTUATIONS IN CASE II 

The quantum fluctuations at initial time are related with the initial data through 



[gj:ni0)r + {qi + a)[gT.um' - = 



TT 
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For > (^^'^^ the positive solution of this equation takes the form, 



where 



1 fgNo^^/^ 



21/3 V 



TT 



2 [gNo I 3 i 



± 



1/3 



For < (^'^^ the positive solution can be written as, 



(2 cos/5-1) 



with 



^qo + (yJ V 27r 



In limiting cases we recover eg. ( |4.14| ). 
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FIG. 1. Case 1. m% > 0. 



36 



Unbroken symmetry, gN^ = 250, qo = 5, Co = 
Comparison between numerical solutions and the early time approximation (|4.7D . 



0, g 



10- 
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qo = 1.3083, a = 0.05233, Co = 0. Comparison between numerical solutions and the early time 
approximation ( ^.18| ). For late times fl'S(r) tends to 1, thus A4%^ = [see Fig. (^)]. 
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FIG. 3. Case II. m| < 0. ^^(t) as a function of r. Dynamically broken symmetry, gNo = 5.101, 
= 1.5336, a = 0.2191, Co = 0- Comparison between numerical solutions and the early time 



approximation 



For late times (not shown in the figure) fl'S(r) tends to 1. 
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FIG. 5. Case II. m'j^ < 0. A^^(r) as a function of r. The symmetry is dynamically broken. 
gNQ = 4.478, qo = 1.3083, a = 0.05233, Co = 0- (The same initial conditions and g as in Fig. §.) 
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FIG. 6. Case I. m|j > 0. Unbroken symmetry. grC{r,T = 2) for g = 10~^ and initial conditions: 
Co = 0, gNo = 250, qo = 4.0, a = 0.3. 
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